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Abstract 

Perturbative QCD relates the single transverse-spin asymmetries (SSAs) for hard processes 
at large transverse-momentum of produced particle to partonic matrix elements that de- 
scribe interference between scattering from a coherent quark-gluon pair and from a single 
quark, generated through twist-3 quark-gluon correlations inside a hadron. When the 
coherent gluon is soft at the gluonic poles, its coupling to partonic subprocess can be sys- 
tematically disentangled, so that the relevant interfering amplitude can be derived entirely 
from the Born diagrams for the scattering from a single quark. We establish a new formula 
that represents the exact rules to derive the SSA due to soft-gluon poles from the knowledge 
of the twist- 2 cross section formula for unpolarized processes. This single master formula is 
applicable to a range of processes like Drell-Yan and direct-photon production, and semi- 
inclusive deep inelastic scattering, and is also useful to manifest the gauge invariance of the 
results. 
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There have been two systematic frameworks to study single spin asymmetries (SSAs) 
observed in a variety of high energy semi-inclusive reactions. One is based on the so-called 
"T-odd" distribution and fragmentation functions with parton's intrinsic transverse mo- 
mentum [1, 2]. It describes SSAs in the region of the small transverse momentum px of 
the final hadrons as a leading twist effect. Its factorization property and the universal- 
ity of "T-odd" functions have been extensively studied [3]. The other one is the twist-3 
mechanism based on the collinear factorization, and is suited to describe SSAs in the large 
Pt region [4, 5, 6]. It relates the SSA to certain quark-gluon correlation functions in the 
hadrons. Recently it has been shown that these two mechanisms give the identical SSA in 
the intermediate pt region for Drell-Yan and semi- inclusive deep inelastic scatterings, i.e., 
describe the same physics in QCD [7, 8]. 

In the present work, we shall focus on the twist-3 mechanism. A systematic study on this 
mechanism was first performed in [5] for the direct photon production p^p 7X, and the 
method has been applied to many other processes such as pion production in pp-collision, 
p^p — > ttX [9, 10, 11, 12], hyperon polarization pp — » A^X [13, 14], Drell-Yan lepton-pair 
production p^p — > £ + £~X [7], pion production in semi- inclusive deep inelastic scattering 
(SIDIS), ep T -> eirX [6, 8, 15]. In our recent work [6], we reexamined the formalism for the 
twist-3 mechanism and gave a proof for the factorization property and the gauge invariance 
of the corresponding single-spin-dependent cross sections, which was missing in the previous 
literature. Through this development, the cross section formula derived in the previous 
studies [5]- [15] have been given a solid theoretical basis. 

In the twist-3 mechanism, the strong interaction phase necessary for SSA is provided by 
the partonic hard scattering: Owing to the insertion of a "coherent gluon" emanating from 
the twist-3 quark-gluon correlation inside e.g. the polarized nucleon, an internal propagator 
of the partonic hard part can be on-shell, and its imaginary part (pole contribution) can 
give rise to the interfering phase that leads eventually to the real single-spin-dependent 
cross section. Depending on the resulting value of the coherent-gluon's momentum fraction 
at such poles, those poles are classified as soft-gluon-pole (SGP), soft-fermion-pole (SFP) 
and the hard-pole (HP). In [6], we have given a formula which expresses all these pole 
contributions in terms of the twist-3 distributions associated with the gluon's field strength 
tensor (see (1), (3) below). There we have also given a general proof that the SGP contri- 
bution appears as both "derivative" and "non-derivative" terms of the twist-3 correlation 
functions, while the SFP and HP contributions appear only as the non-derivative terms. 
Since the derivative of the "SGP function", a twist-3 correlation function at the SGP, causes 
enhancement compared with the non-derivative term in a certain kinematic region, and also 
one may expect that the soft gluons are more ample in the hadrons, phenomenology based 
on the SGP contribution only has often been performed in the literature [9]-[15]. 

Another peculiar feature of the SGP contribution is that the partonic hard cross section 
associated with the "derivative term", arising from the twist-3 soft-gluon mechanism, is 
directly proportional to some twist- 2 partonic cross section as noticed for direct photon 
production [5] , Drell-Yan pair production [7], SIDIS [15], p^p — > nX [9, 10, 11] and pp — > 
A T X[13, 14]. 

In this Letter we propose a new systematic approach to treat the SGPs, which allows 
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us to reveal the novel structure behind the soft-gluon mechanism to the SSA. We disentan- 
gle the coupling of the coherent gluon and the associated pole structure from the partonic 
subprocess, by reorganizing the relevant diagrams using Ward identities and certain decom- 
position identities for the interacting parton propagators. We find that the many Feynman 
diagrams can be eventually united into certain derivative of the Born diagrams without 
the coherent-gluon insertion, which shows that the entire contributions from the SGPs, 
not only the derivative term but also the non- derivative term, can be derived from the 
knowledge of the twist-2 cross section formula for the unpolarized process. We establish 
the corresponding "master formula" that is applicable to a range of processes like SIDIS, 
Drell-Yan and direct-photon production. 

To illustrate our approach, we consider the SIDIS, e(£) +p(p, S±) — > e(t) + ir(P h ) + X, 
following the convention of our recent work [6]: We use the kinematic variables, S ep = (p + 
£) 2 , q — £ — £', Q 2 = —q 2 , x b j = Q 2 / (2p ■ q), and Zf — p ■ Ph/p • q- All momenta p, Ph, I, and 
£' of the particles in the initial and the final states can be regarded as light-like in the twist-3 
accuracy, p 2 = P 2 = £ 2 = £' 2 = 0. As usual, we define another light-like vector as n M = 
(q^ + XbjP^/p-q, and the projector onto the transverse direction as = g^ u —p^n v —p u n^. 
We also define a space-like vector, gf = g M — (P h ■ q/p ■ Ph)p fx — (p • q/p ■ Ph)Ph, which is 
orthogonal to both p and Ph, and its magnitude as qx = \J 1 —q 2 - Then, in a frame where 
the 3-momenta q and p of the virtual photon and the initial nucleon are collinear along the 
z-axis, like the so-called hadron frame [16], the magnitude of the transverse momentum of 

the pion is given by \J —Pyi± = z/9t- 

In the present study we are interested in the SSA for the large- Phj_ pion production, 
in particular, the contribution from the twist-3 quark-gluon correlation functions for the 
transversely polarized nucleon, which are defined as [5, 17] 

= ^ W« ^ PnS± G F (x 1 ,x 2 ) + (7yO y Sj_Gf(x±, x 2 ) + -.., (1) 

where the spinor indices % and j associated with the quark field are shown explicitly, 
F® p is the gluon field strength tensor, and the spin vector for the transversely polar- 
ized nucleon satisfies = g^S^p, Sj_ = —1. The path-ordered gauge-link, \jin, Xn] = 
Pexp [ig fx dtn ■ A(tn)], guarantees gauge invariance of the nonlocal lightcone operator. 
The second line of (1) defines two dimensionless functions Gp{x\ ) x 2 ) and Gp{xi, x 2 ) through 
the Lorentz decomposition of the matrix element; here is the nucleon mass representing 
typical mass scale generated by nonperturbative effects, and "• • •" denotes Lorentz struc- 
tures of twist higher than three. These two functions Gf(xi,x 2 ) and Gf(xi, x 2 ) constitute 
a complete set of the twist-3 quark-gluon correlation functions for the transversely polarized 
nucleon [15, 17]. 

The relevant contributions to the hadronic tensor W^ u arise from the process where the 
partons from the nucleon in the initial state undergoes the hard interaction with the virtual 
photon, followed by the fragmentation into the final state with n + anything, as illustrated 
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(a) (b) 

Figure 1: Generic diagrams for the hadronic tensor of ep^ —>■ eirX, decomposed into the 
three blobs as nucleon matrix element (lower), pion matrix element (upper), and partonic 
hard scattering by the virtual photon (middle). The first two terms, (a) and (b), in the 
expansion by the number of partons connecting the middle and lower blobs are relevant to 
the twist-3 effect induced by the nucleon. 



in Fig. 1. The twist-3 distribution functions contribute to ep 1 ' — > eirX in combination with 
the twist-2 fragmentation function for the pion, which is immediately factorized from the 
hadronic tensor as 



/dz 
—Dj{z)w^{p,q, 
z 



Ph 



(2) 



j=Q,g 



where Dj(z) (j = q, g) is the quark and gluon fragmentation functions for the pion, with z 
being the momentum fraction. We consider the case for the quark fragmentation in detail 
and omit the index j from w^ v below. Modifications necessary for the gluon-fragmentation 
case will be discussed later. The lower blobs in Fig. 1 can be written as Fourier transform 
of the correlation functions for the nucleon, i.e., schematically, M^°'(k) ~ (pS±\%jj%jj\pS±) 
and M^ a (ki, k 2 ) ~ (pS±\'4>A' T ip\pS±) , with the upper indices (0) and (1) representing the 
number of gluon lines connecting the middle and lower blobs. (See Eqs. (23), (24) of [6] for 
the explicit definitions.) The momenta of the partons, k and ki t2 , are assigned as in Fig. 1. 

In the leading-order perturbation theory for the partonic hard scattering, we have 
shown [6] that the twist-3 contribution to w^ v arises solely from Fig. 1 (b), and the en- 
tire twist-3 contribution can be written as 



P h . 



dx\ / dx 2 Ti 



M^(x 1 ,x 2 ] 



dSW(h,k 2 ,P h /z) 



dk% ± 



(3) 



where Mp(xi,x 2 ) is defined in (1), and S^(ki,k 2 ,Ph/z) = Sjp(ki, k 2 , Ph/ z)p a represents 
the middle blob in Fig. 1 (b), which denotes the hard scattering function for the partons 
off the virtual photon, corresponding to the nucleon matrix element (k\, k 2 ). For 

simplicity, we suppress momentum q and the Lorentz indices /i, v for the virtual pho- 
ton in S^(ki,k 2 ,Ph/z). In (3), S^(ki,k 2 ,Ph/z) and Mp(xi,x 2 ) are matrices in spinor 
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Figure 2: Feynman diagrams which give rise to the SGP contributions in the quark frag- 
mentation channel, where the hard quark fragments into the final-state with pion and the 
hard gluon goes into unobserved final state. The cross x denotes the quark propagator 
which gives the SGP contribution. Mirror diagrams also contribute. 



space, and Tr[- • •] indicates the trace over Dirac-spinor indices while the color trace is 
implicit. Here we recall from [6] that, although the straightforward collinear expansion 
to the twist-3 accuracy produces many other complicated terms, it has been proved that 
the hard parts associated with those terms vanish in the leading order in QCD pertur- 
bation theory, using Ward identities for color gauge invariance. Note that (3) guaran- 
tees the factorization property for the twist-3 single-spin-dependent cross section in man- 
ifestly gauge-invariant form, which was assumed in the previous literature [8, 15] (see 
also [5, 7]). We should also mention that, for the HP and SFP contributions, one can 
calculate the partonic hard scattering function via S£\xip, x 2 p, Ph/ z ) by using the re- 
lation dS^(ki,k 2 , Ph/ z)/dk% ± = S^{xip,X2P,Ph/z)/(xi — x 2 ) obtained from the 

Ward identity, while, for the SGP contribution, one has to calculate the "derivative", 
dSW(ki, k 2 , P h /z)/dk% ± , as shown in (3) [6]. 

Based on (3), our task is to identify the SGP contribution to S^(fa, fa, Ph/ z ) an d 
compute its deviation arising linearly in the quark's transverse momentum k 2 ± from the 
value in the collinear limit k± t2 — > x\$p. For this purpose, we work in the Feynman gauge 
and with fa± <C Xtp in fcj = Xip + (ki ■ p)n + ki± (i = 1, 2). In the leading order in QCD 
perturbation theory, S^(ki,fa, Ph/z) stands for a set of cut Feynman diagrams which 
are obtained by attaching the additional gluon to the 2-to-2 partonic Born subprocess, 
where the large transverse momentum Ph±/z of the "fragmenting quark" is provided by 
the recoil from the emission of a hard gluon into the final state. When the coherent 
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Figure 3: The definition of the photon-quark-gluon vertex function T a (k 2) Ph/z), as the 
sum of two types of diagrams which appear in the RHS of the cut in Fig. 2. The amputated 
lines are identified by bars at their ends. 

gluon couples to an on-shell parton line, the parton propagator adjacent to the coher- 
ent gluon produces a pole for the vanishing gluon momentum, k 2 — ki —> 0. Only those 
arising from the diagrams in Fig. 2, where the coherent gluon couples to the final-state 
quark line fragmenting into n + anything, survive as the SGP contributions [15, 6, 8], while 
the other pole contributions cancel out combined with those from the "mirror" diagrams. 
We denote the contributions of the diagrams in Fig. 2, where the coherent gluon is at- 
tached to the LHS of the cut, as S iyl ^ L (k 1 , k 2 , Ph/ z), and those of the mirror diagrams 
as S^ R (k 1 , k 2 , Ph/z), so that S^(h,k 2 , P h /z) = S^ L (k u fc 2 , P h / z) + S^ R (k u k 2: P h / z). 
Explicit form of S <yV)L (ki, k 2 , Ph/z) is given as 

sf> l (k k 2 , ^) = -^f a (k 2 , -)—^f T -—j^—— -iy*!, ^ + h - k 2 ) 

z 2N C z z JDl + ^ - ft 2 + ie z 

xVf(k 2 + q-^), (4) 

z 

where T a (k 2 , Ph/ z) = 7°r^(A;2, Ph/ 2)7° denotes the photon-quark-gluon vertex function of 
Fig. 3 which appears in the RHS of the cut in the diagrams of Fig. 2. In T a (k 2 , Ph/ z), 
the factors for the external lines are amputated. With this T a , the photon-quark-gluon 
vertex function appearing in the LHS of the cut in Fig. 2 is given by r / g(/c 1 , P h /z + ki — k 2 ). 
The structure fh/z projects the final-state quark onto the twist-2 fragmentation process 
to produce the pion (see (2)), and i^aV" is the quark-coherent-gluon vertex. T>°f(k) = 
2n5 (k 2 ) 9(k°)J2\ e f\)(k) e *(\)(k) is the cut propagator for the final-state hard gluon with 
the polarization vector e" x ^(k). Note that we have already worked out the color structure 
associated with (3), simplifying the color matrices contained in Fig. 2 as t a t b t a = (—l/2N c )t b 
and performing the color trace to obtain the gauge-invariant matrix element (1). We also 
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Z 



Figure 4: Diagrammatic representation of Ward identity for the coupling of the scalar- 
polarized gluon to the final-state quark line. 

note the relation S^ R (h, k 2 , P h /z) = S^ L {k 2 , h, P h /z) = -f°S^ L \k 2 , k x , P h /z)^°. 

To proceed further, one may employ direct evaluation of each diagram of Fig. 2, substi- 
tuting the explicit form corresponding to Fig. 3 into the vertex functions T a (k 2 , Ph/z) and 
Tp(ki, P h /z + ki — k 2 ) of (4), and working out the necessary Dirac algebra. Such calculation 
has been done in [6, 8, 15]. Here we employ another approach to disentangle the coherent- 
gluon vertex and the corresponding SGP structure from the diagrams of Fig. 2. An earlier 
attempt in the same spirit investigated decoupling property of the soft-gluon vertex in the 
collinear limit, ki± = k 2 ± = 0, for the SGP contribution in hadron-hadron scattering, using 
the helicity basis technique [20]; but fully consistent evaluation of the SGP contribution 
following [6] requires to treat the deviation of the hard part due to nonzero k 2 ±, as noted 
above, and in this case straightforward application of the helicity basis technique would be 
useless. Our new systematic approach uses Ward identities and decomposition identities for 
the parton propagators interacting with the coherent gluon, which allows us to disentangle 
the coherent gluon vertex from the derivative dS^(k± : k 2 , P h /z)/dk 2± _ in (3). A key 

idea is to reorganize the terms that contribute to this derivative by rewriting p a contracted 
with the quark-gluon vertex in (4) as 

P a = " (*£ - ) " (k a 2± - fcfj.) , (5) 

X 2 — X\ — IE X 2 — X\ — IE 

up to the irrelevant 0((k 2 ± — &i_l) 2 ) correction. Correspondingly, S^ L (ki, k 2 , Ph/z) can be 
decomposed as 

S^ L (k u k 2 , ^) = S^ L (k u k 2 , ^) + S^ L (k u k 2 , ^) , (6) 

z z z 

where the first and second terms in the RHS correspond to those in (5), respectively. In 
(5), "— ie" in the denominator is chosen such that each term in (6) does not produce pinch 
singularity at x\ — x 2 . 

We can exploit an elementary Ward identity in S^ L (k 1 , k 2) Ph/z) in order to disentangle 
the scalar-polarized gluon vertex, ljt 2 — k/i, as well as the quark propagator adjacent to it, as 
(A/*)(#2-MV(A/*+#i-&+ie)] = (P h /z)^ 2 -^-P h /z)[l/(P h /z+^-^ 2 +ie)} = 
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—Jfh/z, as illustrated in Fig. 4, and get 

S^ L (h, k 2 , ^) = -± -f a (k 2 , ^ + h- k 2 )vf(k 2 + (z - — ) . 

z 2N C x 2 — x 1 — is z z z z 

(7) 

This exhibits the SGP only as the single pole, so that we can put x 2 = x\ except for the 
factor l/(x 2 — X\ — ie), without affecting the results for the SGP contribution. Combining 
this with the corresponding result for S^ R (ki, k 2 , P h / z) = S^ L (k 2 ,ki, Ph/ ' z) and Taylor 
expanding the total result S^\k u k 2 , P h /z)= S^ L (k u k 2 , P h / z) +S^ R {k 1 ,k 2) P h / z) in 
terms of ki±_ and k 2 ±, we get 

S^i^hA) = k ^~ K \ ^T a (x lP ,r)^r,(x lP ,r)Vf(x lP + q -r) 

Z 2N C X 2 — Xi — IE OTj_ z 

z 

(8) 

up to the irrelevant terms of the second or higher order in k±±,k 2 ±. Here r denotes a 
four-vector not restricted to being light-like, r 2 ^ 0. 

Next we consider the contribution from S^ 1± ^ L (ki, k 2 , Ph/ z) of the second term of (6). 
Because the coherent gluon vertex associated with S < ^ 1±S)L (k 1 , k 2 , P h / z) is proportional to 
k%± — kf ± (see (4), (5)), we can set ki± = k 2 ±_ = in S^ 1±S)L (k 1 , k 2 , P h / z) except for this 
vertex factor, up to irrelevant corrections. One thus obtains 

S^ L (k u k 2 , = -± T a (x 2 p, ^) 

z 2N C x 2 — Xi — is z 

X —lay: : — -Tp(x lP , — + x lP - x 2 p)Vf(x 2 p + q - — ) . (9) 

Z ^ + (Xi — X 2 )f + IE Z Z 

In this result we decompose the product of the quark-gluon vertex and the quark propagator, 
adjacent to the fragmentation insertion Jfh/z, as 

P h f h jz + {x X -X 2 )i P h ( Pha 1 , Zi \ nn , 

— V2 = — ~d _ + 7<r^5 > l lu ) 

z + ( Xl _ X2 ) p } +l£ z V P h - P x 2 -x 1 -iE 2P h -p) 

where the first term in the parenthesis is given by the eikonal vertex and eikonal propagator, 
and the second term is the "contact term". Combining the result obtained from (9), (10) 
with the corresponding result for S^ 1± ^ R (k 1: k 2: P h / z) = S^ 1±S)L (k 2 , k 1: P h j z), we find that 
the double pole term in S ( ^ 1 - L ^ )L (ki, k 2) Ph/z), which originates from the eikonal propagator in 
(10) and gives the contribution proportional to 8'{x\— x 2 ), cancels the corresponding double- 
pole term in S^ 1± ^ R (ki, k 2 , Ph/z), and the remaining single-pole contributions eventually 
give, for S^KhM.Ph/z) = S^ L (k u k 2 , P h /z) + S^ R (k 1: k 2 , P h /z), 

S^\h,k 2 ,P h /z) = -L ^ - ^, ^ig. P V a {^-)^^-) 

2N C x 2 — xi — ie ar p I z z 



xV^(x 1 p + q- — ) - ^ aPp T a (x lP , r)TfYp(x lP , r)Vf(x lP + q - r) 
Z -i/i ■ p 



(11) 
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Combining (8) and (11), the entire SGP contribution for (3) reads 1 



dSm(h,k 2 ,P h /; 



1 



dS^(x lP ,^) 



ki=Xip 



2N C C F x 2 — X\ — is 



dPi 



(12) 



where C F = {N 2 C - 1)/{2N C ) and 

P h 



,-,/, - . = c F r a ( XlP , ^)^r p(x lP , ^)vf(x lP + q -^) 

z z z z 



(13) 



is exactly the leading-order hard scattering function in the collinear limit k — > x P , which is 
represented as the middle blob in Fig. 1 (a). In (12), we have used the relation, 



9<rp 



PhaP P \ d<p(r) 



Ph-Pj dr p 



M p hl z ) 
dP? 



(14) 



for P h = —(P h 2 ± /2P h ■ P ) P + (P h ■ P )n + P h± , which holds for an arbitrary function <p(r). 
Note that Ph ■ n is not an independent variable to perform the derivative in the RHS of 
(14), corresponding to that the LHS vanishes when contracted by P a . Substituting (1) and 
(12) into (3), and integrating over x 2 to get the imaginary phase from the SGP of (12), one 
gets for w^y as 



1 



w 



2N C C F 



-ttMmz 



apnS± 



d 



dxG F (x, x)Tr 



S^ixpAw 

z 



+ 



(15) 



where the ellipses stand for the contribution from the second term in the RHS of (1), which 
does not contribute to the cross section when contracted with the symmetric leptonic tensor 
for the unpolarized electron. We recall the similar formula for the leading-order twist-2 
hadronic tensor, 

.tw-2 1 f 7 c / \m \ry(n)/ ^ 



w 



[ dxf q (x)Tr S^°\xp, —)]/> 
J L z 



+ 



(16) 



where S^°\xp, Ph/ z) in (15) appears and f q (x) is the unpolarized quark distribution. 

It is straightforward to extend the above results to the gluon fragmentation case by 
mostly trivial substitutions, noting that essential part of our above derivation is based on 
diagrammatic manipulation. The diagrams for the SGP contribution in this case are shown 
in Fig. 5 [6, 15]. The corresponding hard scattering function is given by 



S (1)L 



P h 



iN, 



P, 



Ph 



fei, fe, — = ^K(k 2 , -^)S + (k 2 + g - -±)A p (k u + k 2 ) [-gr(Ph)} 



(17) 



I ki—xi% 



for the SGP contribution, which has 



'We also obtain {d/dk p 1± + d/dk^ ± )S^(k u k 2 , P h /z)\ 
been proved in [6] by the detailed inspection of the diagrams in Fig. 2. This holds for the gluon fragmentation 
channel, too (see (19)). 
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Figure 5: Same as Fig. 2, but for the gluon fragmentation channel where the hard gluon 
fragments into the final-state with pion and the hard quark goes into unobserved final state. 

where the relevant photon-quark-gluon vertex functions A a (k 2 , Ph/ ' z ) and kp{ki, P h / z + 
ki — k 2 ) can be expressed using that of Fig. 3 by appropriate substitutions of the momenta, 
A a (k 2 ,P h /z) = T a (k 2 ,k 2 + q - P h /z). V^^^iqi, q 2 , q 3 ) = (q 1 - q 2 )^g^ 2 + (cyclic per- 
mutation) denotes the ordinary three-gluon vertex except for the color structure that has 
been used as f cab t a t b = (iN c /2)t c to obtain the prefactor iN c /2, the structure —g^ v (Ph) = 
—g ar i + (P^p 71 + PhP a ) I (Ph • p) projects the final-state gluon onto the twist-2 fragmentation 
process to produce the pion (see (2)), and S + (k) = 2n8 (k 2 ) 6{k°)l/ is the cut propagator 
for the final-state hard quark. The mirror diagrams of Fig. 5 gives S^ R (ki, k 2 , P h / z) = 
S^ L (k 2 , ki, Ph/z). Similarly to (6), S^ L can be decomposed as 

S<p L (fc, k 2 , ^) = SW (hi, k 2 , ^) + S™ L (fcx, k 2 , ^) , (18) 

by using (5). The three-gluon vertex coupling to the scalar-polarized coherent gluon in 
S^ L (ki,k 2 , Ph/z) can be disentangled using Ward identity similarly to (7), but this time 
we also obtain some additional ghost-like gauge terms. It is indeed not difficult to show that 
those gauge terms drop in the final result applying Ward identities further. Alternatively, 
one may employ background field gauge [21]: In the Feynman gauge under the background 
coherent-gluon field, the three-gluon vertex in (17) is replaced as V ar] /s(k 2 — k 1 , —P h /z, Ph/z+ 
h - k 2 ) -> V£§(k 2 - h, -P h /z, P h /z + h- k 2 ) = V avf3 (k 2 - h, -Ph/z, Ph/z + h- k 2 ) - 
(Ph/z) v g af3 - (P h /z + ki - k 2 )pg ar) , and the resulting S^ L (k±, k 2 , Ph/z) obeys simple Ward 
identity without unwanted gauge terms, which is analogous to Fig. 4 for the quark-gluon 
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vertex. As a result, we get for S^(k u k 2 , P h /z) = S^ L (k u k 2 , P h /z) + S^ R (k u k 2 , P h /z) 

/ Pi. \ N h a — h a r) 

(huh,—) = -f 2± — ^K(x lP ,r)S + (x lP + q-r)A p (x lP ,r)gf(P h ) 

' (19) 

similarly to (8), up to the irrelevant higher order terms. 

The contribution S^ 1± ^ L (k 1: k 2 , Ph/z) in (18) can be also reduced similarly to (9): We can 
use the decomposition of the product of the three-gluon vertex and the gluon propagator, 
adjacent to the fragmentation insertion gfiPh), m to the eikonal propagator term and the 
contact term as 

9t"( p h)V*p(x2p - x lP , ^ + x lP - x 2P )— -j 

z z (ff + (xi - x 2 ) P ) + ie 

= 9t \P h ) f-^M 1 _ + g avP ") , (20) 

V P h -P x 2 -x l -ie Ph-Pj 

for a =_L. This simple formula analogous to (10) holds in the background field gauge men- 
tioned above. Again, only the single pole terms eventually survive in the total contribution, 
S^iK^Ph/z) = S^ L (k u k 2 ,P h /z) + S^ R (k u k 2 ,P h /z), and the result is given by 
the RHS of (11) with the substitutions necessary for "translating" (8) into (19). Combining 
this result with (19) and taking the similar steps as those in (12)-(15), we find that for 
the gluon fragmentation channel is given by (15) with the replacement 1/N C — > —N c and 
S (0 \x P ,P h /z) -> Sf\x P ,P h /z), where 



Sf\xp, ^) = C F A a (xp, ^)S + (xp + q - ^)Mx P: ^)[-gf(P h )\ . (21) 

z z z z 



This Sf\x P ,P h / z) is the hard scattering function for the 2-to-2 partonic Born subprocess 
leading to the gluon fragmentation, and participates in the twist-2 contribution for the 
unpolarized SIDIS as (16) with S^(x P , Ph/z) — > S^ \x P , Ph/z). 

Substituting (13)-(16) and (21) into (2), we contract the result with the leptonic tensor 
for the unpolarized electron, {£,£') = 2{£ tJi £' ll + £„£') — g^ v Q 2 - Working out the phase 
space factor corresponding to the differential elements, [dui] = dxbjdQ 2 dzfdq^d(f>, with the 
kinematical variables introduced above (1) and the azimuthal angle for the observed final- 
state pion (see [6, 15]), we immediately obtain the formula for the single-spin-dependent 
cross section in SIDIS, e P ^ — > enX, associated with the soft-gluon mechanism induced by 
the twist-3 effects inside the nucleon, as 

da %% _ kM n s± v rdz rdx dH jq {x,z,g*) 

~\M ~ ^ h 9 ' ~' ~ d(P^/z) Gp[X ' X) ' {22) 

where the sum over all quark and antiquark flavors q = u,u,d,d,--- is implicit for the 
index q, and G q F (x,x) denotes the "soft-gluon-pole function" from (1) for the flavor q. 
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The color factors are introduced as C q = —1/(2N C ) and 1/(2N C ) for quark and antiquark 
flavors, respectively, 2 and C g = N c /2. Hj q (x, z,q?) for j = q and g are, respectively, 
equal to Tr S^°\xp, P^/ z)tf and Tr S^(xp,Ph/z)^ , up to the kinematical factor. They 
are exactly the partonic hard scattering cross sections which participate in the twist-2 
factorization formula of the unpolarized cross section for ep — > enX as 

dofw-T 1 ^ r dz r dx 



£ HI -^ D ^)H jq (x,z,q 2 T )f q (x) . (23) 



L 1 j=q,9 

Our results (22) and (23) represent the SGP contribution as the response of 2-to-2 partonic 
Born subprocess to the change of the transverse momentum carried by the "fragmenting 
parton". It is worth noting that our results (22) and (23) hold in any Lorentz frame with 
p± = 0, as is seen from the above derivation. In particular, we recall that the derivative 
with respect to the transverse momentum in (22) was introduced merely as a formal recipe 
via the relation (14), so that one can freely move to any frame even with = after 
performing the derivative. 

In the hadron frame mentioned above (1), P% is parameterized as \J 

Phi. ~ Z flTi 

P^ = ZfQ/y/2 and P^ = —P% ± /2P h ~ = Zfq^/y/2Q, and thus the derivative on the Lorentz- 
scalar functions Hj q (x, z, q 2 -) with respect to P% ± can be performed through qx that is 
indicated explicitly in their argument. Therefore the results (22) and (23) can be expressed 



as 3 



^t S wi _ kM n SiPhi d daZ P 2 l 



(24) 



/,(s)-G« (x,x), Djiz^CjzDjiz) 



[duj] CfZj dq\ [duj 

The partonic Born cross section in (23) was derived as [19, 16, 18] 

= ^rt+rt' (S " (I ~ (i " 0) • (25) 

where Ai — 1 + cosh 2 ip, A2 = — 2, As = — cos0sinh2-?/>, and A4 = cos20sinh 2 -?/>, with 
cosh-?/' = 2xbjS ep /Q 2 — 1, parameterize different azimuthal dependence, and we introduced 
auxiliary variables x = x^/x and z — Zf/z. B 3 k 9 are the functions of x, z, q\, and Q 2 , whose 
explicit form can be found in Eqs. (57) and (59) of [6]. Therefore, the derivative djdq\ 
of (24) can act on either a J k q or the delta-function 5 {q^/Q 2 — (1/x — l)(l/z — 1)), and the 
latter contribution produces the "derivative term" proportional to dG q F (x,x)/dx as well as 
the "non-derivative term" with G q F (x,x), after the partial integration with respect to x. 
We get 3 

ggtgf = a 2 em®se 2 q tiM n s±Ph± 4 v 

dx h ,dQHz s dq 2 d<\> S-Kx 2 bj Sl p Q 2 C F ZfQ 2 ^ k ^ g 3 



2 The relative minus sign in C q between quark and antiquark is due to that of the color charge between 
them, to which the coherent gluon couples. (See Erratum in pp. 18-20.) 
3 See Erratum in pp. 18-20. 
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dG q F (x,x) [l^ 2 <9af 



dx z dq\ 



1 — z dx 



x d(xa 3 k q ) 




(26) 



I 



Substituting the explicit formulae for d k q , it is straightforward to see that this result com- 
pletely coincides with that obtained recently in [6] by direct evaluation of each Feynman 
diagram in Figs. 2 and 5, for all azimuthal dependence, where A\ t 2 give the same azimuthal 
dependence as the Sivers effect [1] while A3 and A4 give additional terms beyond the Sivers 
effect. Our result reveals that not only the derivative term of the SGP cross section but 
the whole partonic hard scattering functions for the SGP contributions are completely 
determined from a k q . 

We emphasize that our result (24), and its immediate consequence (26), show extremely 
nontrivial structure behind the SGP contribution, which would not be unveiled without 
using our new approach discussed above. Indeed it is completely hopeless to infer from 
the complicated formulae (82) and (87) of [6] for the non-derivative term as a function of 
x, z, q\, and Q 2 that the non-derivative term can be expressed in terms of d k q as in (26) with 
a form common to all channel (j — q, g; k — 1, • • • , 4). We also note an important point 
which is most clearly indicated by the compact form (24): Apparently the cross section 
(25) for the 2-to-2 partonic Born subprocess is gauge-independent, and so is its derivative 
with respect to a kinematical variable qr- Combined with (24), this fact guarantees the 
gauge invariance of the hard-scattering function for the twist-3 SGP contribution. In this 
connection we recall that a straightforward check of gauge invariance is hampered for the 
SGP contribution, in contrast to the case for the SFP and HP contributions, because Ward 
identities are useless to the coupling of the zero-momentum coherent gluon. 

The above results (22) and (23) that lead to (24), (26) have been derived using only 
the properties satisfied by the partonic subprocess in QCD perturbation theory. Therefore, 
by analytically continuing the external momenta, the same formula represents the exact 
relations associated with other hard processes. In fact, the corresponding formula for the 
Drell-Yan process can be obtained by switching the final-state "fragmenting parton" into 
the initial-state parton, and the initial-state photon with space-like momentum into the 
final-state photon with time-like momentum. This crossing transformation, as illustrated in 
Fig. 6, is accomplished by the corresponding substitutions: Ph — > —p', 1/z — > x', D q (z) — > 



fq(x'), Dg(z) — > f g (x f ), and g M — > — q^ in (22) and (23), where fq(x') and f g (x') denote 



the twist-2 parton distributions for the unpolarized initial hadron with momentum p', and 
the new q^ denotes the momentum of the virtual photon which is produced by the hard 
interaction of the partons and goes into the lepton pair £ + £~ with the invariant mass 
squared, q 2 = Q 2 , in the final state. With this replacement, our "master formula" (22) also 
describes the single-spin-dependent cross section for the Drell-Yan process p^p — > £ + £~X, 
relating it to (23) for the spin-averaged case pp — > £ + £~ X. The partonic hard scattering 
cross sections are now expressed by the variables, s = (p+p') 2 , s = (xp+x'p') 2 , t = (xp—q) 2 , 



13 



Figure 6: Diagrammatic representation of the crossing transformation of SIDIS into Drell- 
Yan process. 



and u = (x'p' — q) 2 , as (see [7]) 



H jq (s,t,u) 



— — a jq (s,t,u)5[s+t + u-Q 



3ttN c sQ 



(27) 



with (T R = 1/2) 



2C f 



u t 

- + ~ 

t u 



2Q 2 s 
ut 



a° Y (s,i,u) 



2T> 



R 



—is st 



Thus the derivative in (22), which is now with respect to —x'p'^_ instead of P%j_/z, 4 can 
be performed through that for u (see the discussion below (23)). After performing the 
derivative, one can go over to a frame where the two colliding nucleons are collinear along 
the z-axis, and the produced virtual photon has large transverse momentum q±, which is 
provided by the recoil from the hard parton going into the unobserved final state. We thus 
obtains, similarly to (24) and (26), 



da. 



SGP.DY 
tw-3 



a 2 em a s e 2 q 7iM N ^ pnSiQi 



^^/^/^(s + i + u-Q 2 )^') 



dQ 2 dyd 2 q ± 3ttN c sQ 2 C f 



0=1,9 



X 



a dG^pix^ x) 

' -X ; h 



-u 



dx 



a DY 

39 
U 



8a DY 

31 

du 



da 



DY 
39 



t - Q 2 da 



DY-\ 

39 



u ds 



U 



dt 



G q F (x,x)\ ,(29) 



where y is the rapidity of the virtual photon, and Cg = C g . This obeys exactly the same 
pattern as (26) for both derivative and non-derivative terms. Substituting (28), this formula 



4 This corresponds to the fact that the SGP in Drell-Yan comes from the initial-state interaction, while 
that in SIDIS is from the final-state interaction. 
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completely coincides with the result of [7] which was obtained by direct evaluation of the 
Feynman diagrams. 

We can also derive the single-spin-dependent cross section for the direct-photon pro- 
duction, p^p — > 7X, immediately: We make the formal replacement a em /(37rQ 2 ) — > S(Q 2 ) 
in (29), corresponding to the real photon in the final state, and define a® p as the Q 2 — > 
limit of af q Y of (28). Because af q p are invariant under the scale transformation of the 
partonic variables s, i and u, we have (ud/du + sd/ds + id/di)df q p = 0. Accordingly, the 
single-spin-dependent cross section for the direct-photon (DP) production reads 

d 6 q N r s C f J x' J x v ' 



j=q,9 



x d -K { x dG<1 F^ x ) _ G q F (x, x)) , (30) 



with E q = \q\. Note the same twist-2 unpolarized hard cross section appears both for 
derivative and non-derivative terms, as the coefficient for the combination, xdG 9 F (x, x)/dx — 
G 9 F (x,x). This remarkably compact result does not agree with that of [5] 5 , but is reminis- 
cent of the recent result for p^p — > nX [12], where a similar compact formula was obtained 
for the SGP cross section. One can extend the present approach to other processes such 
as p^p — > iiX and show that the compact result found in [12] is also a consequence of the 
simplification due to the scale invariance of the 2-to-2 Born subprocess among massless 
partons, similarly to the present case (30) [22]. 

To summarize, we have studied the SGP contribution in the twist-3 mechanism for the 
SSA. We have developed a new approach that allows systematic reduction of the coupling of 
the soft coherent-gluon and the associated pole contribution, using Ward identities and de- 
composition identities for the interacting parton propagator, and derived the master formula 
which gives the twist-3 SGP contributions to the SSA entirely in terms of the knowledge of 
the twist-2 factorization formula for the unpolarized cross section. We find that this novel 
result is also useful for establishing the gauge invariance of the hard-scattering function 
for the SGP contribution. So far this master formula was derived for DY process, direct 7 
production and SIDIS. Since our diagrammatic manipulation technique uses only elemen- 
tary identities which hold for any diagram with the gluon insertion to cause the SGP, it 
is applicable to other processes such as p^p — > nX and pp — > A^X etc, and may lead to a 
similar relation between the twist-3 SGP contribution and the twist-2 cross section, which 
allows us to reveal the corresponding gauge-invarinant structure. 
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Erratum to: "Master formula for twist-3 soft-gluon-pole 

mechanism to single transverse-spin asymmetry" 
[Phys. Lett. B 646 (2007) 232 (arXiv:hep-ph/0612117)] 

Yuji Koike a , Kazuhiro Tanaka b 

a Department of Physics, Niigata University, Japan 
c Department of Physics, Juntendo University, Japan 



There was an error in Eqs. (24) and (26) of our Letter [1]. These equations were 
derived by performing the derivative in Eq. (22) with respect to the transverse component 
of the final-pion momentum, P% ± = (P^j_, P^ ± , P/uJ = (0, 0,P^j_). The derivative for the 
variation of the magnitude of Ph± was performed through q T in the hadron frame, as 
explained above Eq. (24), but that for the variation of the direction of Ph± also proves to 
yield nonzero contribution. The latter can be performed through associated with A3, A4 
in Eq. (25), where is defined as the azimuthal angle of the initial-lepton's momentum 
measured from the axis along P h± in the hadron frame (see Refs. [6,15] in [1]). The correct 
equations are 



dcr SGP 

uo tw-3 

[duo] 



and 



irM 



N 



C F z) 



c pnS ± P hl _ 



d_ 



-(P 



h± 



dar po1 



'tw-2 



1 d 



2<?T ^0 / [^1 



(24) 



j SGP 
uo tw-3 



liM 



N 



dxbjdQ 2 dzfdq^d(p 8TTxljS 2 p Q 2 C F ZfQ 2 



3=1,9 



c pnS±P h 



fe=i 



X 



dz r dx 



/?/ 



x 



x 



^31 

-jr k x 



dG q F (x,. 
dx 



ji 



dq\ 



x d{xa 3 k 
— z dx 



G q F (x,x) 



x5 



(26) 



+ Q \p 



5, 



dz 

z 



dx 



x 



Dj(z) 



A 2f + A T~ 



G q F (x,x)S\ 




Here the new term with the derivative with respect to in the former formula yields the last 
line in the latter, with As = — sin sinh 2-0 and Ag = sin 20 sinh 2 -0, which parameterize the 
novel azimuthal dependence (see Refs. [16,18] in [1]). We note e pnS±Ph± = Zfqrsm&s and 
(Ph± • S±) = —z/qx cos $s, with $5 that denotes the azimuthal angle for the initial-proton's 
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transverse spin vector S± measured from the axis along Ph± in the hadron frame. From the 
above result (26), it is clear that the twist-3 single-spin-dependent cross section for SIDIS 
consists of five independent structures having different dependences on the two azimuthal 
angles, and $5. We remind that the similar five structures also appear in the SIDIS 
single-spin-dependent cross section obtained by using the transverse-momentum-dependent 
(TMD) parton distribution functions as shown in [2]. Thus, with the new terms in (26), 
the two approaches, our twist-3 mechanism and the TMD one, yield consistent azimuthal 
components for SSA. 

It is instructive to consider the integration of the above formula (26) over the transverse 
momentum of the outgoing pion, P/u_. For this purpose we introduce the azimuthal angles 
4>h and 05 for the hadron plane and the proton's transverse spin vector, respectively, as 
measured from the lepton plane. They are connected to our and $5 as 0^ = and 
05 — — *&s- Then the five different azimuthal structures in (26) can be recast into those 
proportional to sin(0/!— 0^), sin(0/j+05), sin(20^— (ps), sin(30/ l — (ps) and sin 05. Integrating 
(26) over the outgoing-pion's angle 4>h with 05 kept fixed, the contributions associated with 
the first four azimuthal structures vanish, leaving only the term proportional to sin 05. 
However, further integration over q\ (= \P hl _\ 2 / zf) kills this term as well, as can be easily 
checked using the above form (24); the ^-integrals of two terms in (24) cancel out. This 
result is consistent with the fact that there is no SSA in the totally inclusive DIS. 

The results for the Drell-Yan and direct-photon production in the Letter [1], Eqs. (29) 
and (30), which were also derived from Eq. (22), remain unchanged, because those results 
were not concerned with the angular dependence relevant to the above discussion. 

All equations involving the quark-gluon correlation function G q F (x,x) for the flavor q 
in the Letter were presented assuming that G q F (x,y) for the antiquark flavors was given 
by Eq. (1) with the quark fields replaced by their charge-conjugate fields. However, this 
definition appears not to be completely consistent with charge-conjugation invariance. The 
appropriate definition of G q F (x,y) for the antiquark flavors is given by Eq. (1) with the 
quark-gluon three-body nonlocal operator replaced by its charge-conjugate nonlocal oper- 
ator, and in this case all equations involving G q F (x,x), as well as the above formulae (24) 
and (26), should read with C q = —1/(2N C ) for both quark and antiquark flavors. 

The conclusions made in the Letter are unaffected by all the above corrections. 
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